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Renormalization of an Abelian Tensor Group Field Theory: 

Solution at Leading Order 
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Abstract 

We study a just-renormalizable tensorial group field theory of rank six with quartic melonic 
interactions and Abelian group 17(1). We introduce the formalism of the intermediate field, 
which allows a precise characterization of the leading order Feynman graphs. We define the 
renormalization of the model, compute its (perturbative) renormalization group flow and 
write its expansion in terms of effective couplings. We then establish closed equations for the 
two point and four point functions at leading (melonic) order. Using the effective expansion 
and its uniform exponential bounds we prove that these equations admit a unique solution 
at small renormalized coupling. 
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1 Introduction 

Tensor group field theory (hereafter TGFT) is a background-independent formalism for quantum 
gravity. Using the powerful quantum field theory language, it offers both a tentative definition of 
the fundamental degrees of freedom of quantum spacetime and a precise encoding of their quantum 
dynamics. It combines the results of tensor models mta about the combinatorics of random 
discrete spaces and the insights of loop quantum gravity j3j about quantum geometry. More 
in detail, TGFTs are quantum field theories on Lie groups, characterized by a peculiar non-local 
pairing of field arguments in their interactions, whose immediate consequence is that their Feynman 
diagrams are dual to cellular complexes rather than simple graphs. The quantum dynamics is thus 
defined, in perturbation theory, by a sum over such cellular complexes (interpreted as discrete 
spacetimes) weighted by model-dependent amplitudes, in turn functions of group-theoretic data. 
Historically, group field theories (GFTs) [H [5j grew out of tensor models for 3d and 4d gravity 
[6], themselves a generalization of the matrix model definition of 2d Riemannian quantum gravity 
[Tj. In tensor models, the dynamics of a quantum spacetime is given by a sum over equilateral 
d-dimensional triangulations, generated as the Feynman expansion of the partition function for a 
finite rank-d tensor, and weighted by (the equilateral restriction of) the Regge action for simplicial 
gravity. They are thus prototypical models of purely combinatorial random geometries. GFTs arise 
when the domain of the tensors is extended to a group manifold, and the first models [Ij make use 
of these additional data to define amplitudes corresponding to state sum models of topological BF 
theory (by incorporating appropriate gauge invariance conditions, to which we will return in the 
following). Soon it was realized [8j that these group-theoretic data gave the boundary states of the 
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same models the structure of loop quantum gravity states [3]. Later [9], indeed, GFTs were shown 
to provide a complete definition of the dynamics of the same quantum states as their Feynman 
amplitudes are given by spin foam models [lO], a covariant definition of the dynamics of LQG spin 
networks, in turn dual to simplicial gravity path integrals HU- Now, they are understood as a 
natural second quantized formulation of loop quantum gravity [12], and GFT models incorporating 
more quantum geometric features of LQG states and simplicial geometry are indeed among the 
most interesting ones. 

In the meantime, tensor models have witnessed an important resurgence, in the form of colored 
tensor models III El- These solved many issues raised by earlier tensor models and allowed a 
wealth of important mathematical results to be obtained. They triangulate pseudo-manifolds with 
only local singularities [H], having in particular no tadfaces (i.e. a face which runs several times 
through a single edge). Most importantly, they admit a large N expansion [15] (where N d is 
the size of the tensor), whose leading order is now well understood. The leading graphs in this 
limit, the melonic graphs, form particularly simple “stacked” triangulations of the sphere in any 
dimension [TB]. Their appearance is a very general phenomenon mm- Some of these results 
have immediately been extended to topological GFTs and multiorientable models [12], and beyond 
the leading order, to define interesting double scaling limits m- 

Incorporating the insights of colored tensor models into GFTs leads to TGFTs. Here, the 
GFT fields are required to transform as proper tensors under unitary transformations and their 
interactions are required to have the additional U(N)® d invariance, which can be interpreted as a 
new notion of locality, hence singles out a new theory space j'21| . In turn, this invariance requires 
their arguments to be labeled (ordered). Both facts are crucial for GFT renormalization. 

GFT renormalization is in fact a thriving area of current research. Given that the first definition 
of the GFT quantum dynamics is in terms of a perturbative expansion around the Fock vacuum, 
the first aim is to prove renormalizability of specific models, showing therefore their consistency 
as quantum theories. Second, one is interested in unraveling the phase space of GFT models, 
looking in particular for a phase in which approximate smooth geometric physics (governed by 
some possibly modified version of General Relativity) emerges from the collective behavior of their 
pre-geometric degrees of freedom ra. maybe through a process of condensation. The search for 
such a geometric phase, and the associated phase transition(s), is common to tensor models [16] . 
loop quantum gravity [23], and spin foam models [23], but also to other related approaches like 
(causal) dynamical triangulations [23]. Moreover, it has been conjectured [22] to have a direct 
physical interpretation in a cosmological context [26], and some recent results in GFT support this 
conjecture m 

The TGFT framework is well-suited for renormalization, as one can import more or less stan¬ 
dard QFT techniques even in such background independent context. One ingredient is the new 
notion of locality provided by the U(N)® d invariance of tensor interactions. The other ingredient, 
a notion of scale is naturally assumed to be given by the decomposition of GFT fields in group 
representation. This is fully justified in terms of spectra of the kinetic operator (as in standard 
QFT) when a Laplacian on the group manifold is used, as suggested by the analysis of radiative 
corrections to topological GFT models f2S] (which correspond to ultra-local truncations of truly 
propagating models). All these ingredients, it turns out, speak to one another very nicely, as in¬ 
deed in TGFT models counter-terms necessary to cure divergences remain of the same form of the 
initial interactions. More precisely, by precise power counting of divergences, one sees that at large 
ultraviolet (UV) scales (in the sense of large eigenvalues of the group Laplacian) connected sub¬ 
graphs which require renormalization seem local (as defined by tensor invariance) when observed 
at lower scales. 
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A large amount of results has been already obtained. For models without gauge invariance 
the proof of renormalizability at all orders, which started with [29], now includes a preliminary 
classification of renormalizable models [30| and studies of the equations they satisfy [31]. Then 
Abelian [32l 133] and non-Abelian gauge invariance (whose important role we already emphasized) 
has been included PH 35]. The computations of beta functions typically shows UV asymptotic 
freedom mm to be a rather generic feature of TGFTs, even if the analysis of more involved 
models is in fact quite subtle [38]. Renormalizability and UV asymptotic freedom are the two 
key properties of non-Abelian gauge theories which form the backbone of the quantization of all 
physical interactions except gravity, hence it is encouraging to find them also in TGFTs, which 
aim at quantizing gravity. 

Once renormalizability (and possibly asymptotic freedom) is established, the next stage is 
to understand the infrared (IR) behavior of the renormalization group flow, in particular phase 
diagrams and phase transitions. One can prove that the leading “melonic” order of tensor models 
and of topological GFTs exhibits a phase transition, corresponding to a singularity of the free 
energy for a certain value of the coupling mm- The critical susceptibility can be computed at 
least for simple tensor models to be equal to 1/2. In the same tensor models context, in which the 
only notion of distance is the graph distance, one sees a phase corresponding to branched polymers, 
with Hausdorff dimension 2 and spectral dimension 4/3 [4Dj, as in CDT. In GFTs and TGFTs, 
where the group theoretic data play a prominent role, not only computing observables and critical 
exponents, but also finding the nature of the transitions and their physical interpretation is much 
more difficult. 

Therefore we need more analytic tools. One powerful scheme is provided by functional renor¬ 
malization techniques. These have been developed for TGFTs for the first time in [4T]. Applied to 
the (comparatively) easy case of an Abelian rank-3 model, the RG flow equations could be derived 
and the phase diagram be plotted in the key UV and IR regimes, showing evidence for a phase 
transition to a condensed phase, at least in some approximation. 

In this paper we perform a leading order analysis of the correlation functions of a simple TGFT 
with quartic melonic interactions and U(l) group, in dimension 6, endowed with gauge invariance 
conditions. This model is just-renormalizable [33], and asymptotically free Hence it should 
exist at the level of constructive field theory |42| (see [43] for the construction of a simpler super- 
renormalizable TGFT). Although we shall not achieve such a complete non-perturbative analysis 
in this paper, we provide some significant steps in this direction. We define the intermediate field 
formalism for our model and with a multi-scale analysis we establish its renormalizability, compute 
the beta function of the model and check its asymptotic freedom. In this way we recover all the 
results of [33] and [37]. The development of the intermediate field method for our model is in 
itself, we believe, an interesting result. It is known to be particularly convenient for quartic tensor 
models [ 441145] . and should become a standard tool for TGFT’s as well. One should notice in 
particular that in our case, due to the gauge conditions, the intermediate fields are of a vector 
rather than matrix type, a promising new feature. 

We then define the effective expansion of the model, which sits “in between” the bare and 
the renormalized expansion. Its main advantage is to be free of renormalons [42] . We check this 
fact again in our model by establishing uniform exponential upper bounds on effective amplitudes. 
We also establish closed equations for the leading order (i.e. melonic approximation) to the two- 
point and four-point functions. Combining all these results proves that these closed equations 
admit a unique solution for small enough renormalized coupling, and gives full control over the 
melonic approximation of the theory, bringing it to the level of analysis of the Grosse-Wulkenhaar 
non-commutative field theory [45], 
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Similar closed equations have been written for another renormalizable TGFT theory, in dimen¬ 
sion 5 and with a simpler propagator without gauge invariance conditions in [17]. The renormal¬ 
ization and numerical analysis of these equations have been recently developed in |48j. 

Our paper is organized as follows. In Section 2 we define the model and its intermediate field 
representation. In Section 3 we establish and analyse its power-counting with multi-scale analysis. 
Section 4 describes its renormalization, computes the beta function (in agreement with 133). 
introduces the effective expansion and establishes uniform bounds on the corresponding effective 
amplitudes. Section 5 writes the closed equations for the melonic approximation to the bare and 
renormalized two point and four-point functions, and completes the proof that these equations 
have a unique solution at small renormalized coupling, which is in fact the Borel sum of their 
renormalized expansion. 


2 The Model 

In this section, we shall briefly recall the basics of TGFTs models with closure constraint (gauge 
invariance) and Laplacian propagator. Then we shall focus on a particular U( 1) quartic model 
at rank six first defined in [32] • Within this section definitions and computations are still formal 
since we do not introduce cutoffs; this will be done in the next sections. 


2.1 General Formalism for TGFTs 


A generic TGFT is a statistical held theory for a tensorial held, for which the entries are living in 
a Lie group G, generally compact, such as 1/(1) or SU 2) for the simplest cases. A family of such 
models was defined and renormalized to all orders in [32], [34, [33] I35Q 

The theory is defined by an action and by the following partition function 

S(0,0)=S int (0,0)~J-0~0-J, Z{J,J) = J dn c $, (2.1) 


where Si nt is the interaction and d/i<y is a Gaussian measure characterized by its covariance C. The 
helds 0 and 0 are complex functions 0, 0 : G d C noted 0(gi, • • • ,9d) — 0(g), g = (gi, # 2 , gd), 
and 0(g'), g 1 = (g],^, ■■■,9d)- They should equivalently be also considered as rank-d tensors, that 
is elements of the tensor space L 2 (G)® d , where L 2 (G) is the space of functions on G which are 
square-integrable with respect to the Haar measure. The 2A-point Green functions are obtained 
by deriving N times with respect to sources J and N times with respect to anti-sources J 


G 2 n{9 


i, 


i 9ni 9 1,''' 9 n) ~ 


d 2N z(J, j ) 


dJ^g^dJiig’J ■ ■ ■ dJ N (g N )dJ N (g'j 


N) 


J=J =0 


( 2 . 2 ) 


The Gaussian measure is defined by the choice of the action’s kinetic term. TGFTs such as 
those of [22] '3t)J use a mass term plus the canonical Laplace-Beltrami operator A on the group G d , 
hence correspond to the formal normalized measure 


d/dcM 0 ) = -^-e~ Skin D(f>D(j) (2.3) 

4 Renormalizability has not been yet established for models based on the Lorentz group, which is non-compact. 
However, at least intuitively, one could expect the additional difficulties present in the non-compact case to be 
rather of IR nature than of UV nature, from the point of view of TGFT renormalizability; this would imply similar 
renormalizability results as in the compact group case. 
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with 


S t in{ 0 , 4 >) = / [dg] d <p(g)[(-A + m 2 )^](g), 


(2.4) 


where dg is the Haar measure on the group. Although the Lebesgue measure D(f>D(f> in (2.3) is 
ill-defined, the measure dpc 0 itself is well-defined, and the propagator Co in the parametric (or 
Schwinger) representation is 


Cu(g,g")= / oj«®S! , rj = / Aae a ’ n ‘ fj I< a (g c g' c '), 


(2.5) 


C=1 


where K a is the heat kernel associated to the Laplacian operator, and c is our generic notation for 
a color index running from 1 to d. In momentum space this propagator becomes diagonal. Let us 
from now on restrict to the case G = 1/(1). The Fourier dual of 1/(1) is Z, hence in momentum 
space, we note p = (p 1? • • • ,Pd) G Z d , where p c G Z is called the strand momentum of color c, and 
we have 

d i 

c 0 ( p , p ') =nwj „, ( 2 - 6 ) 

C—1 


p 2 + m 2 


In the specific TGFT we study in this paper, we want the field configurations to obey the 
additional gauge invariance 

0(Pi, 92 , ■■■, gd) = 4>(hgi, hg 2l ..., hg d ), 0(pi, g 2 , ...,g d ) = 4>(hg 1 , hg 2l ..., hp d ) V/i G G. (2.7) 


This gauge invariance complicates slightly the writing of the model. In order to implement it, we 
could introduce the (idempotent) projector P which projects the fields on the subspace of gauge- 
invariant fields, then equip the interaction vertices and propagators with such projectors. But 
in this case the tensorial symmetry U(N)® D symmetry of the interaction vertex (which provides 
the analog of a locality principle for renormalization) would be blurred. Hence the best solution, 
used in [32], consists in implementing the gauge invariance directly on the Gaussian measure by 
introducing a group-averaged covariance 

/ /•OO f d 

dpc(0,0)0(p)0(p') = y dae~ am2 J dhYlKaigchg'c 1 ). (2.8) 


In other words, we introduce the gauge invariance projector P only in the propagator of the 
theorjj^] In momentum space we have. 


d 

C\P,P') = Yl^PcPc) 

C—i 


HJ2 c Pc) 

p 2 + m 2 


(2.9) 


From now on we shall remember that the covariance is diagonal in momentum space, with diagonal 
values 

C(p) = 6 ^_ cP j - (2.10) 

pZ _|_ m z 

hence defining the set V = {p G Z 6 | J2 c Pc = 0} of momenta satisfying the gauge constraint, 
all Green functions of our theory can in fact be defined for restricted momenta p G V, or if one 
prefers, are zero outside V. 

5 Additional insertions of P on the vertex would result in the same Feynman amplitudes, since P 2 = P. 
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TGFT interactions by definition belong to the tensor theory space mmm spanned by 
U(N)® d invariants. Hence the most general polynomial interaction is a sum over a finite set B of 
such invariants b , also called d-bubbles, associated with different coupling constants tb 

S in t(4>, <t>) = hh{4>, <t>), (2.n) 

fees 

where If, is the connected invariant labeled by the bubble b. Graphically, each bubble is associated 
with a bipartite rl-regular edge-colored graph. Each color c G {1,2,.., cZ} is associated with a 
half-line at each vertex, and each vertex bears respectively a field (j) or its complex conjugate (f) 
according to its black or white color. The edge coloring of the bipartite graph allows to visualize 
the U(N)® d invariance by showing the exact pairing of fields and anti-fields argument of the same 
color. Such graphs also enable to visualize whether the interaction is connected or not. Some 
examples of connected invariants at ranks d = 3 and d = 6 are shown in Figure [lj 



Figure 1: Some connected tensor invariants 


The Feynman amplitudes of the perturbative expansion are associated with Feynman graphs 
whose vertices belong to the set B of the interaction d-bubbles. A Wick contraction is represented 
by a dotted line. Figure [2] gives an explicit example for d = 3. 



Figure 2: A tensorial vacuum (N=0) rank-three Feynman graph 


For a Feynman graph Q, we note V(£7), £{G) and £{G) the sets of the vertices (the d-bubbles), 
internal (dotted) lines and external (dotted) half-lines, and V(Q), L(Q) and E(Q) = 2N(Q) the 
number of elements in these sets. The number of vertices V is also identified with the order of 
perturbation, also often noted n. 

The Green functions are given by a sum over Feynman graphs (connected or not) 

G 2 n= WHY AG' ( 2 ’ 12 ) 

g,E(G)= 2 N \beB J 
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where rib is the number of vertices of type b and s(Q) is the graph symmetry factor (dimension 
of the automorphism group). Note that expanding each vertex b as a d-regular bipartite edge- 
colored graph as in Figure [T| and coloring the dotted lines with a new color 0, any such graph Q is 
therefore canonically associated to a unique (d+ l)-regular bipartite edge-colored graph, for which 
the vertices are the black and white nodes, as shown in Figure [2] Hence it defines an associated 
d-complex, in which in particular faces are easily defined as the bi-colored connected components 
BUSS!. These faces are either closed or open if they end up on external half-lines. 

The connected Green functions or cumulants G'^y are obtained by restricting sums such as 


(2.12) to connected graphs Q, and are obtained from the generating functional 


=log[Z(J,J)] 


(2.13) 


through 


G 2 n( 9 i j ‘‘ 5 9 n-> 9 n''' 9 n) ~ 


d 2N w(J , j) 


dJ^dJ^g'f) ■ ■ ■ dJ N (g N )dJ N (g' N ) 


j=j =o 


(2.14) 


The vertex functions r 2j y are obtained by restricting sums such as (2.12) to one particle irreducible 


amputated graphs Q (amputation mean we replace all the external propagators for dotted half-lines 
by 1). They are the coefficients of the Legendre transform of W(J, J). 

Using the convolution properties of the heat kernel (following from the composition properties 
of its random path representation), the Feynman amplitude Ag of Q can be expressed in direct 
space as 


Ag — 



dhe 


n a «.v 

f&HG) 



II K °m 


9*(/)f l Kel h ‘' 9 



(2.15) 


In this expression, F(Q) is the set of internal faces of the graph, F ext {Q) the set of external faces, 
and etf the adjacency matrix which is non zero if and only if the line £ belongs to the face / and 
is ±1 according to their relative orientation. We noted «(/) = the sum of Schwinger 

parameters along the boundaries-lincs of the face /, and g s ^) or g t ^ the boundary variables in 
the open face /, s for “source” and t for “target” variables. We use also the notation T for the 
set of faces and F for its cardinal (number of elements). 

These amplitudes Ag can be interpreted as lattice gauge theories defined on the cellular com¬ 
plexes dual to the Feynman diagrams Q. The group elements hi (resp. g s (f), 9t(f)) define a discrete 
gauge connection associated to the edges £ (resp. boundary edges) of the cellular complex, and 
the ordered products ( res P- 9s(f)Y\. e &df^ 1 9t{f)) are hs holonomies (discrete curvature) 

associated to bulk (resp. boundary) faces of the same complex]^ 

Due to the diagonal character of the propagator in momentum space, these Feynman amplitudes 
are easier to express in the momentum representation. In particular the momentum conservation 

6 In models of 4d quantum gravity that bear a closer relation with loop quantum gravity, and that encode 
more extensively features of simplicial geometry, additional conditions called simplicity constraints are imposed 
m\mm- Obviously, they complicate the structure of the amplitudes, making them richer. We do not consider 
these additional constraints here. 
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along faces due to the 5 functions in (2.9) ensures that when expressed in momentum space non¬ 
zero Green functions of the theory of order 2 N must themselves develop into sums over U(N)® d 
tensor-invariants of the momenta of order N ; in other words to any entering momentum p c must 
correspond an exiting momentum with same value p' c = p c . In particular the two point function 
in momentum space is a function of a single momentum p G Z d , and the connected four 

point function G1 is a sum over all quartic invariants of the theory. In general the contribution 
of a given specific tensor invariant is complicated to extract from the Green functions. It requires 
a somewhat subtle decomposition using Weingarten functions, which we shall not detail here, 
referring the reader to 


2.2 The Quartic Melonic U(l)-model in dimension 6 

After this quick overview of general TGFTs, we come to the particular model studied in this 
paper, namely the d — 6 Abelian quartic model with melonic interactions. It is the simplest 
just-renormalizable model (with no simplicity constraints) in the classification of gauge invariant 
TGFT models [31]. As such, it is also the simplest interesting testing ground for the analytic 
techniques we develop here. 

General quartic interactions at rank 6 are of the three types indicated in Figure [3] Melonic 
interactions correspond to the type 1. They are leading in the 1/IV tensorial expansion and are 
marginal in the renormalization group (RG) sense, the other ones being irrelevant. 




Figure 3: The quartic tensor interactions at rank 6 


Hence the interaction part of the action considered from now on is the sum of all the bubbles 
of type 1. There are 6 of them, characterized by a unique index c referring to the special color 
which colors the two lonely lines of the bubble: 


6 

S int = Yl A cTRc(#)- (2-16) 

C—1 

More explicitly a quartic interaction b\ with special color 1 writes 

Tr bAH) = I dgdg'$(gi,g 2 r" ,96)<f>(g[,92r" ,9s)<j>(9i,92r’- ,9*) 

= ; P6)0(Pi,P2,-- - ,P6)0(p1,P2,--- jPe), (2T7) 

p,p' 


where the last line is written in Fourier space. Remark that since only fields satisfying the prop¬ 
agator constraints Pc = 0 can contribute, in (2.17) we must have p\ = p \. Hence each Tr ftc is a 







function of fields with 9 (rather than 10) independent strand momenta, because p c = p' c . We can 
therefore in our model simplify ( |2.17[ ) into 

Ti'ft, (#) = $(PuP2,--- ,Pe)(i>(pi,P2,--- ,Pe)(i>(pi,P2,--- ■ ■ ■ ,Pe)- ( 2 - 18 ) 

p£'P,p'£'P\pi=p l 1 

From now on we consider only the color-symmetric case A c = A Vc = 1, • • • ,6. 

As remarked, Green functions in momentum space develop into sums of tensor invariants. In 
particular the connected four point function G 4 develops over all quartic invariants (connected 
or not). Hence it develops over the connected invariants of Figure [3] and over the disconnected 
invariant which is the square of the quadratic invariant. This may seem dangerous at first sight 
since to be renormalizable our model should not involve in particular renormalization of invariants 
of type 2 and 3 which are not part of the initial interaction. 

As well known, renormalization is best stated in terms of the vertex functions T. Hence we 
shall be particularly interested in computing the two point vertex function or self-energy Y 2 (p) 
and the four point vertex function F 4 {p\ , ...,p 4 ). These functions are a priori defined on V or V 2 . 
However we shall see that their divergent part is simpler. More precisely we shall define melonic 
parts YVf elo (p) and rf iel °(p l , ...,p 4 ) for these vertex functions, and even a refined monocolor melonic 
part Y™Y ono (p c ,P'c) of YT d {p\ , ...,p 4 ), such that Y 2 (p) -Y™ elo (p) and Y^ mono (p c ,p' c )^Yf^° ono (p c ,p' c ) 
are superficially convergent (hence truly convergent after all divergent strict subgraphs have been 
renormalized). More precisely we shall prove that 

Theorem 1 There exist two (ultraviolet-divergent) functions f and g of a single strand momentum 
p6Z such that 

6 

rr ‘>1 = -A Er (p„p'Mpc)- (2.19) 

C— 1 


and such that Y 2 {p)—Y™ el ° (p) and Y^ mono (p c ,p() — Yf l ff° no ( y p c ,p' c ) are superficially convergent (hence 
truly convergent after all divergent strict subgraphs have been renormalized). All higher order vertex 
functions are also superficially convergent. 


In particular Y 


melo o nr ] T^melo 

dliu 1 4 mono 


(Pci Pc ) both depend in fact of a single non-trivial function, respec¬ 


tively / and g, of a single strand momentum in 7L. We shall prove that the special form (2.19) of 
the primitive divergencies of the theory is compatible with the renormalization of the couplings in 
(2.18). In the next section we introduce the intermediate field representation in which the functions 
/ and g are particularly simple to represent graphically and to compute. 


2.3 The intermediate field formalism 

The intermediate field formalism is a mathematical trick to decompose a quartic interaction in 
terms of a three-body interaction, by introducing an additional field (the intermediate field) in the 
partition function. It is based on the well-known property of Gaussian integrals: 

/+00 

dxe~ x2/2 e iKxy = V^e~ K2y2/2 . (2.20) 

-OO 

We first apply the general method without exploiting gauge invariance, then stress the simplifica¬ 
tion due to gauge invariance. This means we start with (2.17) which we want to exhibit as a square. 
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For this we introduce the six auxiliary matrices ^ 2 P2 ... P6 4>(pi,P2, • ■ ■ ,P&)4 > {p'\iP2i ■ • • ,pe) = M pi p ', 
which are quadratic in terms of the initial 0 and 0 and can be thought as partial traces over color 
indices other than 1. The interaction in (2.17) can be rewritten as 


Tr 6l (00) = tr 


( 2 . 21 ) 


where tr means a simple trace in £ 2 {Z). Using many times (2.20) we can decompose this square 
interaction tr M 2 with a new Hermitian matrix o\ corresponds graphically to “pinching” the two 


special strands of color 1 with this matrix held, as indicated in Figure |4j More precisely 

f daie~ tr ( <T i)/ 2 e *v / 2A tr (oiM) 


—Atr ( 


2 ) _ 


f da 1 e tr Gi)/ 2 


( 2 . 22 ) 



Figure 4: Intermediate held decomposition 


The next step is to make this decomposition systematic for the six melonic interactions. Writing 

tr (<7 iM) = Tr0S 1 0, (2.23) 

where 

S 1 = cri(8)I(8)I(8)I(8)I(8)I (2.24) 

acts in the large tensor space £ 2 (Z) 06 and Tr means a trace in this large tensor space, allows to 
express the previous intermediate held decomposition as 

A tr M 2 J daie~ tr (° r l)/ 2 e i v / 2ATr V 

f d<7ie~ tr Gi)/ 2 

Using color permutation, we decompose all six bubbles in this way. An intermediate held a c is 
therefore associated to each quartic bubble b c with weak color c. The operators 

E c = I-. • 0 a c <g) ■ • -I (2.26) 

commute in the tensor space £ 2 (Z) lg)6 , as they act on different strands. Introducing E = X^c=i ^ 
we can rewrite the partition function of the original theory as 

Z(J,J) = j d/zc(0,0)e- 5 “ t( ^ ) = y r d/i C (0,0)di/((x)e-^ J -^ J V 72XTV ^, (2.27) 

the normalized Gaussian measure dv(a) being factorized over colors with trivial covariance iden¬ 
tity on each independent coefficient (Gaussian unitary ensemble). The tensor integral becomes 
Gaussian, hence can be computed as a determinant. We hnd: 
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Proposition 1 The partition function of the model is given by 


Z(J,J) = / e-^-^XC^CJ-TrHl-iV^C^^y 


(2.28) 


Therefore pairs of sources are become resolvents (1 — iy/2\CTi) 1 in this representation. 
Perturb at ively one can expand both the interaction logarithm and these resolvents as 


- Tr ln(l - iy/2XCT) = -yTr (iV2XCE) n -, (1 - iV2\CT)~ x C = Y^(iV2\CT) n C. (2.29) 


n =1 


n= 1 


We call the factors Tr ( i\/2XCTj) n loop vertices [50] and the factors (i\/2\CY>) n C ciliated vertices 
HU or, more simply, chains. 


We now incorporate the important simplification (2.18) due to the gauge constraint of our 


model. It ensures that all components of the a matrices factorize trivially from the integral (2.28) 


except the diagonal ones. More precisely since for an intermediate matrix of color c only the p c = p' c 
term contribute, any loop vertex or chain depends only of the diagonal part r c (p c ) := (cr c ) PctPc of the 
six intermediate field matrices previously introduced, ffence we can reduce the six intermediate 
matrices in our model to six vector fields r (these diagonal parts)Q Since each r operator is 
diagonal we conclude also that all propagators occurring in either a Tr (z\/2ACS) ri loop or in a 
(■ i\/2XCTj) n C chain have the same momentum p e 7L d . Since we remarked that the E c operators 
all commute together in the tensor space £ 2 (Z)® 6 , the value of a loop vertex is a simple sum over 
the numbers k\, ■ ■ ■ ,k$ of insertions of er 1 , ■ ■ ■ , a 6 , their total number being n. It can therefore be 
written as 


n\ 


s (EcPc 


n 


Tr (iV2\CY,) n = [iV 2A] n V r ^ . 

pez 6 , fceN 6 1 J2 c k c= n ilc=i c c = l 

P&V 

where we recall that C 0 (p) = (p 2, + m 2 ) -1 , V = {p e Z 6 | J2 c Pc = 0} an d we define 

Tip) = V2X ^T c (p c ). 


(2.30) 


(2.31) 


Similarly any chain is a diagonal operator, hence depends on a single momentum p and is non-zero 
only for p G V, with value 


\ C / fceN 6 I y] c fc c =n ilc=1 c c = 1 

= [iC 0 (p)T(p)] n C 0 (p). 


Y[ T c(Pc 


(2.32) 


Hence 

' This important simplification could be interesting for a future constructive analysis of the model. 
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Proposition 2 The partition function of the model is given by 


Z(J , J) 


J dv{r)e~^vev ^(p)( 1 -jC' 0 (p)T(p)) 1 c 0 {p)j{p)-J2 j r er H 1 - iC o{p)T(p)\ _ 


(2.33) 


where dv is the normalized Gaussian measure on the six vector fields r c (p), each defined on 7L, with 
trivial covariance 

J dv(T)T c (p)T d (p') = 5(c,c’)5(p,p'). (2.34) 

We want now to describe graphically the Green’s functions G 2 N and the vertex functions of 
the initial theory in this intermediate field formalism. 


2.4 Graphical Representation 


This subsection provides our graphical conventions and Feynman rules for the intermediate field 
perturbative expansion in the momentum representation. An intermediated field propagator is 
represented by a wavy line, which bears a color label c and carries a single momentum p c , and 
correspond to the covariance of a r c intermediate field. The loop vertices (which come from deriving 


the logarithmic interaction in (2.28)), are represented by grey disks, to which intermediate field 
half-lines are attached. The chains (which come from deriving the source term in (2.28)) are 


represented as ciliated lighter gray disks: they are characterized by a single ciliurn, represented as 
a dotted half-line attached to the disk, see Figure [5} A ciliurn has no color and represents on its 
left side the entrance of the particular momentum p of the chain and on its right side its exit. 



Figure 5: Intermediate held graphic representation: propagator, loop vertex and ciliated vertex 


The former propagators C which were the dotted lines of the initial representation of Figure [2] 
are now in one-to one correspondence with the arc^ on the boundary of all the disks (both the 
loop vertices and the ciliated vertices), see Figure [bj 

Green functions G%n of the initial theory can be computed as Feynman graphs with exactly 
N ciliated vertices and an arbitrary number of loop vertices H- In particular G 2 correspond to 
the sum over connected graphs with exactly one ciliated vertex, and G C A to the sum over connected 
graphs with exactly two ciliated vertices (a generic one is pictured in Figure [b]). 

We can also consider the pure intermediate held theory with the J and J sources put to zero. 
We introduce new sources J dual to r. A source J is therefore a set of six functions J c (q c ) for 

8 These arcs are often called corners in the mathematic literature; here we prefer a more physical terminology to 
convey the fact that arcs are associated to propagators. 
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Figure 6: A 4-point graph of the tensorial theory in the intermediate field representation 


q = {q c } G Z d . Introducing the natural notation J ■ r = Y2qez d ^2c=i ^c(q c ) ■ r c (g c ) this pure 
intermediate field theory is dehned by the partition function 


Z(J) = 


di i (r)e~ J " r ~^-‘p^ v T(p)\ _ 


(2.35) 


It has connected Green functions corresponding to expectation values of products of r fields 


W£'- CN (q u --- ,q N ) = 


d N log z(j) 


dj ci (?i), ■■■dJ CN (q N ) 


J =o 


(2.36) 


These expectation values are represented by a sum of Feynman graphs such as those of Figure [7], 
with a total number of q wavy half-external lines attached to the loop vertices (grey disks), each 
carrying a color c and a single strand momentum q c . 



Figure 7: A 4-point graph of the pure intermediate field theory 


By color permutation symmetry, the one point function of the pure theory at color c, W^(q^) 
is in fact independent of c. We call it therefore W\{q). It is a function on Z. Similarly the pure 
intermediate fields two point function q 2 ), which a priori is given as a function of two 

colors c and c! 1 and of two strand momenta q\ and q- 2 , can by color permutation symmetry be 
described by just two functions on Z 2 , namely W^(q±, q- 2 ), which corresponds to ci = c 2 , and 
Wf(qi,q 2 ), which corresponds to C\ ^ c 2 . 
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The renormalization of our model will involve only the melonic approximation of the 2- and 


4-point vertex functions Tg ° and r™ el '°. But there is a simple correspondence between the melonic 
approximation of the 2iV-point vertex functions T^° of the initial theory and the tree approxima¬ 
tion of the TV-point Green functions W^ ee of the pure theory, discovered in the context of tensor 
models m- We shall develop it in our case in subsection |3.3| and use it to identify graphically 
the functions / and g in Theorem [T} We now return for a while to the initial theory to establish 
its power counting and renormalization using a multi-scale analysis. This analysis will lead us 
naturally to focus on the melonic approximations which govern renormalization. 


3 Regularization and power counting 

3.1 The regularized theory 

Simpler superrenormalizable Abelian TGFT models [32] as well as a just renormalizable non- 
Abelian model at rank 3 [MJ have been analyzed already using a multiscale expansion. We recall 
the basic steps of that analysis here, adapting it to our specific model. 

Like any theory with ultraviolet (UV) divergencies, this model requires a UV cutoff before 
introducing the renormalization procedure (which gives a coherent scheme to extract finite and 
cut-off independent information). We shall use in this paper both the parametric cutoffs as in [35] 
and sharp momentum cutoffs, which are simpler for our model because of its strong momentum 
conservation rules. 

The parametric cutoffs slice the Schwinger parameter. We fix a parameter M > 1 and define 

Co (g,g’) 

Ci(g,g') 

We choose the UV-regulator A so that A = M _2p , and the complete propagator C\ = C p is then 
given by: 

p 

= (3.2) 

i =0 

A corresponding sharp momentum cutoff X<p(p) is 1 if |p| 2 < M 2p and zero otherwise. The 
theory with cutoff p is defined by using the covariance 

C p (p)=C(p)x< P (p)- (3-3) 


'i 


I M~ 


I 

dae ~ am2 / dh\[K Q (g c hg'- 1 ) 

J c= 1 

n d 

d ae ~ am2 / dhJjKaigJig'; 1 )^ ± 0. 




(3.1) 


C—1 


Then we slice the theory according to 

c e ($) = J2Clip), C,(p) = C(p)Xi(\p\ 2 ) (3.4) 

i =1 

where Xi is 1 if \p \ 2 A M 2 and zero otherwise and for i > 2 Xi is 1 if M 2( *W < |p| 2 < M 2,1 and 
zero otherwise. 


14 



A subgraph S C G in an initial Feynman graph is a certain subset of lines (propagators C) plus 
the vertices attached to them; the half-lines attached to the vertices of S (whether external lines of 
G or half-internal lines of G which do not belong to S ) form the external lines of G. Translating to 
the intermediate representation, we find that a subgraph should be a set of arcs of the intermediate 
field representation, plus all the wavy edges attached to these arcs. The external lines are then 
the (half)-arcs attached to these wavy edges which do not belong to S. 

A vertex of the initial representation is called external for S if it is hooked to at least one 
external line for S. Similarly a wavy line of the intermediate representation will be called external 
to S if it hooks to at least one external arc. 

Particularly interesting subgraphs in the intermediate field representation are those for which 
the set of arcs are exactly those of a set S C CV of loop vertices (excluding any chain, so no arc 
belongs to any ciliated vertex). Let us call such subgraphs proper intermediate or PI. Remark that 
any PI graph is automatically 1PI in the initial representation (since all arcs belong to at least 
one loop, the one of their loop vertex). Also any PI graph can be considered amputated, hence 
as a graph for a particular vertex function. The converse is not true and many graphs for vertex 
functions do not correspond to PI graphs in the intermediate representation. 

PI subgraphs can be represented as graphs of the pure intermediate theory, simply by omitting 
the two half-arcs at the end of each external wavy line. In our model their amplitude depends 
only of the single strand momentum entering the wavy line, not of the full momentum of the two 
half-arcs hooked at its end. 

We shall see that in our theory only very particular non-vacuum connected subgraphs are 
superficially divergent, namely PI graphs which are trees with at most two external lines. 

3.2 Multiscale analysis 

The multi-scale analysis [123 allows to renormalize in successive steps, in the Wilsonian spirit. It 
attributes a scale to each line £ E C{Q) of any amplitude of any Feynman graph Q. 

Let us start by establishing multi-scale power counting. We can perform this analysis both 
with parametric or sharp cutoffs, ending with the same conclusions. In this subsection we use the 
sharp cutoffs since they attribute the same scale to all arcs of any loop vertex or chain, hence a 
single scale to any loop vertex of the intermediate field representation. 

The amplitude of a graph Q , A(Q), with fixed external momenta, is thus divided into the sum 
of all the scale attributions p = {ie,£ G £(C/)}, where ie is the scale of the momentum p of line £: 

= < 3 - 5 ) 

At fixed scale attribution p, we can identify the power counting in powers of M. The essential 
role is played by the subgraph Qi formed by the subset of lines of Q with scales higher than i. 
By the momentum conservation rule along any loop vertex, this subgraph is automatically a PI 
subgraph which decomposes into k(i ) connected PI components : Qi = These connected 

components form when (i, k ) take all possible values an abstract tree for the inclusion relation (the 
famous Gallavotti-Nicolo tree [49 ]). We have 

Theorem 2 The amplitude A^(Q) is bounded by: 

\Afj,(Q)\ ^ K l{q ) n fl M^\ K > 0, (3.6) 

i k= 1 
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and the divergence degree io{fH) of a connected subgraph ji is given by: 

u(U) =-2L{H) + F(H) - R{H), (3.7) 

where L(fH ) and F{fH) are respectively the number of lines and internal faces of the subgraph PL, 
and R{fH) is the rank of the adjacency matrix e^ for the lines and faces of PL. 

Proof: Obviously we have (for K = M 2 ) 

\C,m < K6Qjp c )M- 2i x<i{p). (3.8) 

C 

Fixing the external momenta of all external faces the Feynman amplitude (in this momentum 
representation) is bounded by 


\MS)\ < 


[j KM~ 2it 

leC(G) 


n En x <^p> n <*(&*)■ 

/£F,„(6)i>/EZte»/ te£(0) 


(3.9) 


(k) 

The key to multiscale power counting is to attribute the powers of M to the Q\ connected 
components. For this, we note that, trivially: M l = M -1 J([*. =0 T7, a trivial but useful identity 

which allows e.g. to rewrite rifeL(G) M~ 2h = M 2 rifeL(G) ECLo ^ Then, inverting the order of 
the double product leads to 


k(i) k(i) 

J] M~ 2ie =n n M ~ 2 =nn n M~ 2 = nn M~ 2 L( Fi\ (3.10) 

l£L(g) i i£C ( u‘Wgt) * k=1 i&C(g^) i k =1 

The goal is now to optimize the cost of the sum over the momenta p / of the internal faces. Summing 
over pf with a factor x<i(p) leads to a factor KM 1 , hence we should sum with the smallest values 
i(f) of slices i for the lines £ G df along the face /. This is exactly the value at which, starting 
form i large and going down towards i — 0 the face becomes first internal for some Qf. Hence in 
this way we could bound the sums n feF int (g) 


k (‘i) 

YIYI mF(g ' ] - ( 3 - n ) 

i k =1 

However this can be still improved, because we have not yet taken into account the gauge factor 
n ee£(g) 3(I2cP c )- It clearly tells us that some sums over pf do not occur at all. How many 
obviously depends of the rank R of the incidence matrix e^f. Indeed rewriting the delta functions 
in terms of the Pf(e,c) we have 


II 5 E^) = II 5 (J2 e ifPf)- (3T2) 

eec{g) c i&c{g) / 


Hence writing the linear system of L equations "jffj e efPf — 0 corresponding to the delta functions 
we can solve for R momenta pj in terms of L — R others. It means that in the previous argument 
we should pay only for F — R sums over internal face momenta instead of F [^} 

9 Remark that the remaining product unused or redundant S functions are simply bounded by 1 because the pf 
variables are discrete, hence the S function are simply Kronecker symbols, all bounded by 1; of course this would 
not be true for continuous variables as a product of redundant 5 distributions in the continuum is ill-defined. 
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This argument can be made more precise and rigorous and distributed over all scales starting 
from the leaves of the Gallavotti-Nicolo tree (the smallest subgraphs Qf) and progressing towards 
the root we can select faces such that the restricted sub-matrix q/ still has maximal rank R(G-) in 
each Qf. We discard the other faces decay factor. Then we can select lines so as to find a restricted 
square submatrix egf with maximal rank R(Gi) in each Gf■ This leads to 

k(i^ 

| A^G)] <K c ^ g) nn m -2 L(e?)+F(Gj)-fl(Sf) = K \L(G)\ nn (3.13) 

i fc=l i k =1 

This equation completes the proof, and the exponent u)(Gi) = —2L(Gi) + F{Gf) — R(Gf) identifies 
the divergence degree. □ 

3.3 Melonic graphs 

In this subsection, we will determine the nature of PI superficially divergent graphs, which are 
those with positive divergence degree oj > 0. We shall establish that they are melonic |1] in the 
ordinary representation, and trees in the intermediate field representation. 

Consider first the case of a PI vacuum subgraph. If it is a tree on n loop vertices, it has 
L = 2 (n — 1) arcs, 5n + 1 faces (since each wavy line glues two faces) and it is easy to check by 
induction (adding leaves one by one from a root) that the rank R of the e matrix is maximal, 
namely n. Hence oj = —4 (n — 1) + 5n + 1 — n — 5 in this case. 

Next let us consider the case of a PI tree subgraph with N external wavy lines, hence 2N 
external arcs. 

• if N = 1 the subgraph is a two point function and the single external wavy line adds one 
arc, suppresses one face and does not change the rank, hence oj — 2 in this case. 

• if N — 2 the subgraph is a four point function and the two external wavy lines adds two arcs. 
If they have different colors, or have the same color c and hook to two components of the 
tree not connected by lines of color c, then they open two different faces and do not change 
the rank, so that oj — — 1. However there is a special case, when the two external wavy lines 
have same color and hook to the same loop vertex or to different loop vertices joined by a 
path in the tree made of wavy line all of the same color c. In that case and only that case, 
the wavy lines open only the single face of color c common to all loop vertices along this 
path, the rank again has not changed and oj — 0. 

• if N > 2, each new external line takes L into L + 1, can either keep F unchanged (if it hits 
an already open face), in which case R is also unchanged, or takes F to F ■— 1, in which case 
either R is unchanged or goes to R — 1; Hence oj decreases at least by 1. This proves 

oj(G) < —(N - 2) if N > 2. (3.14) 

Consider next the case of a PI vacuum subgraph with N external wavy lines and q wavy loops. 
We can first pick a tree of wavy lines then add the wavy loops one by one. Each added loop creates 
two new arcs and changes the number of faces by -1, 0 or 1. It can change the rank at most by 1, 
and when it creates a face, then the rank cannot decrease (the matrix e becoming bigger). Hence 

oj(G) < -(IV - 2) - 3 q if N > 2. (3.15) 
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In particular if N = 1 and q > 1 we have 00 (G) < —1 and the graph is convergent. 

Finally it remains to study the case of non-vacuum, non-PI graph. Since they add at least 
one new arc to a PI graph, it is easy to check they have lu < 0, except in two particular cases 
corresponding both to one-particle reducible graphs: 

• a chain of arcs joining PI two-point trees, with one of them at both ends. Such subgraphs 
are one-particle reducible two point subgraphs of the initial theory with uj — 2. 

• a chain of arcs joining PI two-point trees, with one of them at a single of its two ends. Such 
subgraphs are one-particle reducible four point subgraphs of the initial theory, with to = 0. 

These cases are not interesting since such subgraphs cannot occur as and, as is well known, 
renormalization can be restricted to IPI subgraphs. 

These results in particular show that the degree of divergence c <j does not depend on the number 
of vertices, but only on the number of external lines. This is typical of a just renormalizable field 
theory. 

Trees in the intermediate representation correspond to melonic subgraphs in the ordinary rep¬ 
resentation JITj. Hence we have proved, in agreement with the other renormalizable TGFT’s: 

Theorem 3 The only superficially divergent PI subgraphs are melonic in the ordinary representa¬ 
tion, with two or four external ordinary lines. In the intermediate representation, amputating the 
trivial external arcs, they are PI trees with a single external wavy line, or with two external wavy 
lines of the same color carrying the same strand momentum. 

Melonic graphs are graphs with zero degre<^_ J hence for which all jackets are planar. We include 
for completeness brief definitions of these two notions, referring to [Ij for more details. 

Definition 1 (Jackets) A jacket J of a regular d + 1 colored graph Q c is the canonical ribbon 
graph associated to Q c and to a (D-hl)-cycle f up to orientation. It has the same number of lines 
and vertices than Q c , but contains only a subset of the faces, those with consecutive colors in the 
cycle Tj = {fe T Qc \f = ififiifii), ^ +1 (0)), q G Z D+1 }. 

Hence there are dl/2 jackets at rank d and to each jacket is associated a Riemann surface of 
genus gj. 

Definition 2 (Degree) The degree w{Qf) is by definition the sum over the genus of all the jack¬ 
ets: 

w (Gc) = =>- m(Gc) > 0- (3.16) 

J 

The degree governs the 1/N tensorial expansion since the number of faces is a monotonically 
decreasing function of the degree. Melonic graphs have maximal number of faces at a given 
perturbation order. More precisely 

Lemma 1 The number of faces F c of Q c is related to the number of black vertices p and to the 
dimension d by: 

did — 1) 2 

F c = - 2 — ~P + d ~ ( d _i)\ w ^c)- (3.17) 

10 The degree in question is the ’’degree of the colored graph”, which characterizes the dominant order of the 
large-N limit of tensor models. It should not be confused with the degree of divergence, and we denote it by w(Q c ). 
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A tensorial graph Q having a unique colored extension Q c , we can extend the notion of degree 
to tensorial graph. Since the colored extensions of type 1 vertices of our theory all have the same 
number of inner faces (faces without color 0), the degree of Q c again governs the number of faces 
of Q, which are the bicolored faces of Q c which includes color 0. In our case the vertices of Q all 


have 25 inner faces and p = 2 black vertices, so that (3.17) tells us 


F(Q) = 5V + 6- -—w[C,c 
60 


(3.18) 


2 and Tf el ° are naturally defined as the melonic approximations to U and 


Returning to Theorem [3] we can precise the divergent part of the theory in the language of the 
previous section. Yf el ° and T™ eio 
r 4 and Theorem [3 indeed proves that U — T 
Moreover they express simply as tree approximations of the pure r intermediate field theory: we 
have 


■^meio an( ] y 


4 ,mono ~ ^Tnmno are superficially convergent. 


rr'°w = v'ZX Yw?“(p c ), TT“(p c ,p' c ) = s/2\S(p c ,p' c )W 2 - ,ree (p e ,p c ) (3.19) 


c= 1 


where W[ ree and Wf’ tiee are respectively the tree approximation to W\ and Wf. 

But Theorem |3] contains still an additional information on the divergent part of rf el °. Defining 
W 2 ~mono as the P ar t °f W 2 ,ree in which all wavy lines along the unique path joining the two external 
lines must be of the same color c than these two external lines, it states that the difference 
is also ultraviolet finite, hence can be neglected in the following section on 


W« 


Wo 


2 vv 2 ,mono 

renormalization. 


Since (3.19) is nothing but (2.19) with / = —i 
the proof of Theorem [lj 


tree.c 


and g — W. 


tree 

2, mono ? 


this completes 


3.4 Uniform Convergent Bounds 

An important aspect of the multiscale analysis is that it provides easily a uniform exponential 
bound on convergent amplitudes: 

Theorem 4 (Uniform Weinberg theorem) The amplitude A(Q) for a completely convergent 
connected graph Q (i.e. a graph for which co(H) < 0 \/TL C Q) is uniformly bounded in terms of its 
size, i.e. there exists a constant K such that if n is the order (number of vertices) of the graph: 

\A(G)\ < K n{G \ (3.20) 


Proof We assume N(Q) > 1, so that MTL C Q,N{fH) > 1 (the vacuum case N(Q) = 0 is an easy 


extension left to the reader). (3.14) implies that for a convergent PI graph with 2N > 4 external 
arcs 

u{H) < -N(H)/3 = -2N(H)/G. (3.21) 

This is also true if TL is convergent with N = 1 or 2, since we saw that in this case cu < — 1 < 
—N(fH)/ 3. For a (j) A graph of order V = n with 2N external legs, we have 2 L = 4V + 2N. Therefore 


(3.13) implies that for another constant K 


k{i) 


A{g) < k u znn M 


-2AT(0f)/6 


(3.22) 


k= 1 
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Let us now define 


L,(/i) = sup ie(n),e v (n) = inf iRyu), (3.23) 

£e l„(S) teL 6 (5) 

where u denotes a vertex v G £/, and L V (Q ) the set of its external (half)-lines. v is external to a 
high subgraph C/f if and only if e& < i < %, and then it is hooked to at least one of the 2 N(Q k ) 
external half-lines of Q k . Therefore 

JjM“ 2JV( ^ fc))/6 < Yl II (3.24) 

vegj k) \e v <i<i v 


Using the fact that there are at most 4 half-lines, and thus 6 = 4 x 3/2 pairs of half-lines hooked 
to a given vertex, and that, for two external lines i and i' of a vertex v, \e v — i v | > | in — iy\, we 
obtain: 

A(G)<K n J2 II n (3-25) 

H v (g/')±v 


where the product over (£,£') _L v means the product over all pairs of half-lines hooked to v. 
The bound means that there is exponential decay in scale differences between all such pair, t ]~ n ~[ 
Organizing the sum over /i = along a tree of lines of Q as in [42], it is easy to bound it by 
, hence to complete the proof of (3.20), hence of Theorem |dj □ 

The next section is devoted to renormalization of the model and to a computation of its beta 
function. 


4 Perturbative Renormalization and Flow 

Renormalization consists, after having identified the “dangerous” subgraphs Q k (those with uj > 0), 
in subtracting from them their local Taylor approximation (the “counter-terms”), up to cancelation 
of the divergencies, hence up to order oj. Then one should compute how renormalization changes 
the interaction from bare to renormalized, hence compute the flow of the theory from the ultraviolet 
to the infrared regime. 


4.1 Perturbative renormalization and counter-terms 

Our goal in this section is to check that, as stated in [35j 37J 

Theorem 5 The U( 1) model with T 4 interaction at rank 6 is just renormalizable and asymptoti¬ 
cally free. 


The perturbative renormalization implies the following redefinitions 


0 = Z x !\ A)0 r , 0 = Z l '\ A)0 r , 


(4.1) 


A = Z-\K)Z]l\K)\ r = Z]! 2 A r , (4.2) 

m = Z -1 / 2 (A) Z]l 2 (A)m r = Zl( 2 m r , (4.3) 


n 


(3.25) is of course a very sloppy estimate, that could be easily improved. For instance we could take advantage 
of the momentum representation conservation rules to remark that only one pair of different scales is in fact hooked 
to any vertex, rather than 6, but it won’t change the structure of the result, only improve numerical constants. 
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and the UV-regularized generating partition function is: 


Z: = J ^ c(z _ 1/2z v 2mr) (Z 1/2 0 r ,Z 1/2 0 r )e^ 2 ^^= 1 ^^^) 

= [ dfi c{z _ 1/2z ij 2rnr)/z {<j> r , f r )e z ^ K TnMrtr)' 


(4.4) 


In these definitions, the “r” subscript applies to the renormalized quantities. The mass and 
wave function counter terms can be absorbed in the covariance 

j d/ J ' C (Z-y 2 zK 2 m r )/Z^ r ^ ^r)0r(^)0r(^) 


with 5z — Z — 1, 5 m 2 = Z m — 1. Identifying this covariance with the one of the initial bare theory 
means that the (bare) propagator of the theory rewrites in terms of renormalized quantities as 




s(E c Pc 


jp-(e-d') 


^ Zp* + Z 


-l 


Z m ml 


E 




P 2 + m 2 r 1 5 z f + S m 2 m 2 r 


0 ip(e-d') 


(4.5) 


j 9 


mi 


C 


d(J2 c Pc) _ 

p 2 + ml + SzP 2 + Supm 2 ' 


(4.6) 


As well known the renormalized parameters in a BPHZ scheme are obtained in terms of the 
bare ones through the vertex functions, which are the one-particle irreducible amputated functions. 
In our model the power counting analysis of the previous section showed that we need only to 
renormalize the 2 point vertex function T 2 (self-energy), and the four point vertex function T 4 . 


4.2 The renormalization group flow 

The basic idea of the renormalization group is the following. All correlation functions are invari¬ 
ant under an infinitesimal dilatation s 1 + 8 of the ultraviolet cut-off A with a simultaneous 


redefinition of the coupling constants, mass, and field normalization: 

A —y A(1 + 5), m —> m + Sm, A —> X + 6X, Z^-Z[1 + 5Z). (4.7) 

Renormalized quantities parametrize a given trajectory of the RG flow. We have the relations: 

4> s a = Z 1 ^ 2 (sA)(j) r , m(sA) = Zll 2 {sA)m r , \{sA) = zlj 2 (sA)\ r , (4.8) 

involving 

fsA = Z 1 / 2 (sA)Z~ 1 / 2 (A)f A =: Z^ 2 (s)<f > A , (4.9) 

m{sA) = zK 2 (sA)Z~ 1/2 (A)m(A) =: Z^ 2 (s)m( A), (4.10) 

A(sA) = zl /2 (sA)Z; 1/2 (A)X(A) =: ^ /2 ( S )A(A). (4.11) 


These relations give the transformations of field, mass and couplings of two theories with dif¬ 
ferent cut-offs, hence along the same trajectory of the RG flow. They imply trivially the invariance 
of the renormalized correlation functions along a given trajectory. This invariance translates into a 
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differential equation for the correlation functions describing the evolution of the RG flow, namely 
the so called Callan-Symanzik (CS) equation. Writing 


G^ mAl ;jm) = z N {K)G^ Xr m\ 


with 


[Z(sA)Z-\A)] N Gl N mAM = G™ msA:XsA . 
and developing, with s := 1 + <5, we get 

Z(sA)Z~ 1 (A) = (z(A)+A^-s) Z~\ A) = l + 5A-^lnZ, 
\ dA J dA 

n 2N — r™ A A j — | ^A d dm" 8 | 2N 

^sA,m sA ,X sA - ^A,m A ,X A + ° dA 8\ dA 8m 2 J 

Gluing the pieces, we obtain the CS equation: 

A 7)A + + m2 ' ym2 ^drr^ + ^(A) j G S» a ,Aa = °> 


dA v '8X 

with the following definitions: 

r 2N ._ }_ 

^A.rn A ,X A ■— ^ 


N 


ducM <t>) n 

3 = 1 


n a \ ^ i r? \ d i 2 

/?:=A dA''' := - A dA lllZ ' 7 '” 2:=A dA lnm ' 


(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 


(4.17) 

(4.18) 


We analyze now this equation at first order (one loop). 


4.3 One loop self energy 

We start by computing the corrections to the propagator. At one loop, the only melonic (hence 
divergent) graph is pictured in Figure [8] 




Figure 8: The self-energy contribution at one-loop 


Its value is 


-pmelo,l 

1 2 


ip) 


EE*a 


( q c ) 2 + m 2 


E 

c q c c , cYc 


(g^) 2 + m 2 ’ 


(4.19) 
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where q° G 7L d has components q c d = (qi, q 2 ,q c — P, • ••, q<i)- 

Let’s use the Schwinger representation to rewrite the denominator as an integral of an expo¬ 
nential, with UV cutoff on the parameter a. ft gives 


/ OO 

- /A5 <' 




- 2A E 


dcie" 


l/A 2 


, 2vr 4^ ' 




q c'^c 


dae~ 


^ d(3 . 


2vr 


e if3pc e -ap; 


nE 

c'A 


e i/3<Z c c , e -a(<£) 2 


(4.20) 


In the last equality, we introduced an integral representation of the Kronecker delta. Now, we can 
turn the argument of the exponent into a perfect square and obtain, for p e V 


rr ,o ’ i (p)=- 2 AE r dae ~° m * rz 

Jl/A 2 Jo 271 


d[3 _kr2 


e i(3p c e ~apl e ~hf3 2 /4a 


nE 

cYc g‘ 


e -a( q ‘,-il3/ 2a) 2 


-2AV [ dae~ am2 [ ^ e i/3pc e- Qp ' e - 5/32/4 " (-) 5/2 , 
^./i/a 2 Jo 2vr W 


(4.21) 

(4.22) 


in the a —>■ 0 limit. This identifies the divergent behavior of this expression. These divergencies 
come from the neighborhood a — 0, and using the distributional expansion 


e -p 2 /4a = + a 5"((3)\ + 0(a 5 /2), 


(4.23) 


we obtain: 


rr^tp) = -2A V [ dae ~ Qm2 [ —e ifipc e- ap2 V^[S(fJ) + ad” (£)] (-) ^ + 0(1/A) 

, Jl/A 2 Jo 2 7T 

= “5 L L A “ r “ v “ Pl (i "*4) + ° (1/A) - (4 - 24) 


The asymptotic expansion of this expression at large A is now easy to find using integrating 
by parts 


Z = 


/■°° 2 2 1 
/ dae~ am e~ ap c—. 

'l/A 2 « 


'l/A 2 


dae~ am --pi / dae ~ am “—b C?(l/A) 

a J i/a 2 a 


= AV m2/A2 - (p 2 c + m 2 ) / dae-“ m2 - + 0(1/A). 

./l/A 2 « 


(4.25) 


The divergent part of the last integral is at most logarithmic near zero. Thus: 


/ dae~ am — = A ln(A) + 0(1/A), 

1 1/A 2 Cl 


(4.26) 
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and it suffices to determine A. Differentiating with respect to A and identifying the singularity in 
the two expressions, we find 


4r r dae~ am2 - = -^e“ m 2 /A 2 A 2 = \e ~ m2/A2 


dA 


' l/A 2 


a A 3 


A 


= | + 0( 1/A 3 ) = A(/3)i + 0( 1/A 2 ) => A = 2, 


and we obtain the following divergent part: 

I = A 2 — 2 (p 2 + m 2 ) In (A) + 0(1). 


(4.27) 


(4.28) 


Returning to (4.24), we find then 


-pmelo, 1 

1 2 


(p)(p) = -^=- J2 ( a2 _ 2 ^ 2 + m ' 2 ) ln ( A ) _ I7 2 ln ^ A )) + °( 1 / A ) 

12A7T _ 2 m 2 in (A)) + ln(A)p e + 0(1/A). 
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5^5 


(4.29) 


and comparing with (4.6) we conclude that at one loop 

24Att 2 


— 


<S m 2 77r = 


575 

12Avr 2 

7H~~ 


hi(A), 


(A 2 — 2m 2 In (A)) 


(4.30) 

(4.31) 


4.4 Coupling constant renormalization and asymptotic freedom 

In this section we examine how the coupling changes along the RG trajectory i.e. going towards 


the IR.. Equations (4.29)-(4.30) gives us the coefficient 7 at first order: 

d 24Av t 2 

7 = -A— lnZ - - 

dA 575 


(4.32) 


imelo, 1 


It remains now to evaluate the melonic monocolor four-point function at one loop, T r 4mono . 
The contributing diagram is sketched in Figure [9] (shown with its four external arcs). 



Figure 9: One loop melonic 4-point function 


The total contribution of this diagram is 


1' 


>£ £ 

c <f c l c'/c 


7E C '7') 

[(g ^) 2 + m 2 ] 2 


dm 2 


a £ £ 

c c 7 c 


7E C ' 7') 

(<f c ) 2 + m 2 ’ 


(4.33) 
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and we can deduce the divergent part of X' (in the same notations as in the previous section): 

T = MA) + 0( 1/A). (4.34) 

The last thing to evaluate are the symmetry factors. We have four ways to connect the external 
fields in an amputated vertex, and two ways to connect these contracted vertices for each of the 
6 colors. The expression of the four point function at zero momentum is then ultimately, to the 
(leading) one-loop order: 


/ 0\2 2 \ 

r Ci = 6 x 4 (-A + -j=- ln(A) j + 0( 1/A). 


Returning to equations (4.16) and (4.18) we get: 


2A 2 tt 2 


/3(A) 


48A 2 7t 2 

5\/5 


0, 


(4.35) 


(4.36) 


which implies immediately, 


a £ = « a > 


38A 2 tt 2 

5^5 


(4.37) 


The minus sign is fundamental. It means that the bare coupling constant decreases when the 
ultraviolet cutoff increases. The theory is therefore asymptotically free, thus consistent at the 
perturbative level, like the familiar non-Abelian gauge theories of the standard model. 

We now discuss what happens beyond one loop. 


4.5 Counterterms and renormalons 

Renormalized amplitudes A^Q) can be explicitly written in terms of Zimmermann’s forest formula 

Ar(S)= Y. (4-38) 

F&D(Q) 7 €-F 

where r* is an operator which performs explicitly the subtraction of the counter-term and D(Q ) is 
the set of all divergent forests of Q. However such renormalized amplitudes suffer from the problem 
of renormalons. Indeed they can grow as n\ with the number n of vertices in Q. This problem 
exists also in our model, and even in its mclonic approximation. Consider indeed the two point 
subgraph of Figure [Tdj made of an arbitrarily large monocolor chain of n simple loop vertices with 
two arcs, ending on a leaf with a single arc. All wavy lines have same color c and carry the same 
momentum p c . 



Figure 10: Typical melonic graph with renormalon effect 


Because the renormalized four point function, hence the renormalized loop vertex with two 
arcs, behaves as log(p c ) at large p c , inserting such a chain on a convergent loop in a convergent 
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melonic vertex function will lead to a very large sum over p c which typically can behave at large 
n as 

~ K " A (4.39) 

p c 6Z 


P 


2 _j_ m 2 


for some constant K. This is the renormalon problem. 

Snch renormalons in fact come entirely from the counterterms in Zimmermann’s formula [ 42] . 


More precisely in (4.38) the counter-terms are subtracted, so to speak, blindly with respect to 
internal scale integrations. But a divergent subgraph looks like its counter term only when its 
internal scales are higher than the scale its external lines, hence when it is a Qf in some attribution 
/i in the language of the previous section (locality principle). Counter-terms in (4.38) not only 


subtract these dangerous contributions, but also the inoffensive parts in which the internal lines 
of the divergent subgraph have lower scale than the external lines. It is exactly these unnecessary 
subtractions which give rise to the renormalons. Hence, although the standard renormalization 
procedure eliminates all ultraviolet divergencies from any Feynman amplitude, such renormalized 
amplitudes are so big that we cannot use them directly to sum even the melonic approximation to 
our theory. 

4.6 The effective amplitudes 

The effective series is a more physical way to compute perturbation theory, and a natural solu¬ 
tion to the renormalon problem when the theory is asymptotically free [42]. The basic idea is 
to renormalize in the Wilsonian spirit, namely step by step, expanding in a whole sequence of 
effective couplings rather than in the single renormalized coupling. Consider a graph Q and its 
bare amplitude A^Q) with scale attribution p as defined in the previous section. There are some 
Qf subgraphs which are divergent (uj(Qf) > 0). They form a forest D^{Q) (because it is a subset 
of the Gallavotti-Nicolb tree containing all G( high subgraphs). The effective amplitude A e ^(G) 
is defined by 

= A l"<.5):= n (1 -r*)A.(S). (4.40) 

v 7eny(S) 


Comparing with (4.38), we see such amplitudes are very different from the renormalized ones. 


Because in (4.40) all divergent high graphs are subtracted, effective amplitudes, like renormalized 
ones, have a finite limit when the ultraviolet cutoff is removed. However unlike renormalized 


amplitudes, effective amplitudes are free of renormalons 


More precisely 


Theorem 6 The effective amplitude A e fj{G) for a graph Q with V(G) internal wavy lines is uni¬ 
formly exponentially bounded in term of its size, hence for some constant K 


\Ae fJ (S)\ < K Vi0> - 


(4,41) 


Proof (sketched) Renormalization operators exactly act on the divergent subgraphs Gf only. Taylor 
expanding and using the condition that external legs of G) have all lower scales than any internal 
line, they transform their divergent degree into an effective degree a/ < —1. The rest of the 
argument to bound the sum over p then follows exactly the proof of Theorem [4} □ 
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Hence effective amplitudes are better building blocks than either bare or renormalized ampli¬ 
tudes to understand the ultraviolet limit of the theory. It remains to relate them to the initial 
theory. Consider the bare power series defined by: 


n ~ 


G,n 


s(Q) 


■X) vm A^9), 


(4.42) 


where attributions /i are summed with cutoff A = M~ 2p (hence by (3.4) every scale satisfies 
1 < i( < p) the amplitudes A^ are computed with bare propagators and A is the bare coupling. It 
has obviously no ultraviolet limit. But we have the following key theorem 


Theorem 7 (Effective expansion) The series (4.42) can be reshuffled as a multi-series with 
effective couplings and effective amplitudes: 


A 

U 2N 


Z 


s(S) 


n (- a 


(A) 

Am) 


iuev(6) 


>j Af(S), 


(4.43) 


where i v (p) = sup {%iff hooked ton}, and the effective couplings and the effective propagators 

C/ ff occurring for lines of scale i in the amplitude AjfffQ) obey the inductive relations (4.45)-(4.46) 
below. Moreover, defining the renormalized coupling by X r := Ao and the renormalized propagators 
C r by inverting (4.6), and reshuffling the effective series in terms of the single renormalized coupling 
X r with renormalized propagators C r , we recover exactly the renormalized series. 

In particular Ap A ' ) is the bare coupling, and Aq A ^ is the renormalized one. The other couplings 
A* for 0 < i < p describe the RG trajectory in between these extremal values. 


Proof (sketched) We recall only the main steps in the proof; more details can be found in [421 S]- 
The proof is inductive, working from the high scales towards the lower ones. The initial step 
i = p starts with the bare series. At step number i we suppose we have defined the effective 
expansion with 

• effective couplings A j for vertices with highest scale j > i and A i+ i for all vertices with highest 
scale j < i] 

• effective propagators Cj- for lines with indices j > i and Cj+i for all lines with indices j < i, 

• effective amplitudes A e ^' l+1 {Q) with subtractions n 7 eD)+ 1 (6)^ — w h ere D^f 1 is the forest 
of all divergent Q 1 ) with j > i. 

We define the next coupling \ and propagator C % by considering in /i the scale number i. 
Adding and subtracting the counter-terms in \ = {'H e DffQ) \ inf i e — i),l e TL, we write 

Af’ i+ \S)^ n [(l-r;) + rj] n (! - (4.44) 

neD^Dff 1 i&Dff 1 

and we expand the product over TL G D 1 \ D l ff l . The operators (1 — rffl) will generate the next 
layer of subtraction in the formula to change the subtraction operations of A e ff^ l+l {Q) into those 
of A e Jf ,l (Q). The counterterms +r^ are then associated to collapsed graphs Q/TL in which TL has 
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been collapsed to a vertex (if N(H ) = 2) or to a mass or a wave function insertion (if N(H) = 1). 
Collecting these pieces and rearranging them according to the collapsed graph rather than to the 
initial graph defines an (infinite series) redefinition of the couplings hooked to vertices with highest 
line of slice j < i and of the propagators with scale j < i, which become respectively A* and C % . 
Hence the new effective coupling is 

-A, = -A j+1 + V T ’ n Af’‘+ >(W) (4.45) 

H\ N('H=2),inie G - H i£=i 

and the new propagator is 

Ci = C,+i+ V T^A e "' <+1 m (4.46) 

H\ N(H=2),m{ e fz n ie=i 


Remark we can omit in these definitions that H is divergent, since = 0 if "H is convergent. 
Remark also that H in (4.46) is connected but can be one particle reducible and that to update 
the effective mass and effective Laplacian normalization in C r from i + 1 to i requires to analyze 
(4.46) in terms of the one-particle irreducible self-energy (see (4.6)). Finally remark also that such 
recursive equations are non-Markovian. By this we mean that the effective coupling A, is itself 
a multi-series in the sequence of all effective couplings A p , - - • , A,;+ 1 , Any attempt to rewrite it in 
terms of the single coupling Aj+i would automatically reintroduce the renormalon problem. □ 
Thanks to Theorem [6] the effective expansion is therefore able to define the theory provided 
all couplings on the trajectory from A p to Ao = A r are uniformly bounded by a sufhcently small 
constant, and the number of graphs is not too big. This is the case when 


• the theory is asymptotically free or asymptotically safe in the ultraviolet regime, 

• the set of graphs considered does not proliferate more than exponentially with size n. 

Planar “wrong sign” 0 4 [5Tj or the Grosse Wulkenhaar model [46] satisfy these two conditions. 
Since melonic graphs, like trees, obviously proliferate no more than exponentially in size and since 
our theory is asymptotically free, its melonic approximation also satisfy both conditions. Hence 
the effective expansion allows to define non perturbatively this melonic approximation, in fact for 
any Green function G^§ 1 ° or vertex function T^°. 


5 Melonic Equations 

In this section, we establish a closed equation for the melonic two-point vertex function, and an 
equation expressing the melonic four-point vertex function in terms of the two-point one. Com¬ 
bining this with the effective bounds of the previous section we shall prove existence and unicity 
of the solution of these equations at small renormalized coupling. 

5.1 Bare Equations 

Let us start with the two-point vertex function or self-energy. The relationship between the Green 
function G 2 and the self-energy can be graphically represented as: 
corresponding to the functional relationship: 

c . S(Z c Pc) 

1 — r 2 c p 2 + m 2 — r 2 (p) 


g 2 (p) = c + cr 2 c + cr 2 cr 2 c + ■■■ 


(5.1) 








g 2 



C , c ( \ C , C C ( c , 

+ -( 1PI - +- -^-4 1PI - +■ 




Figure 11: Decomposition of the 2-point function 


We want to restrict now this relationship to the melonic approximation. (3.19) expressed Y™ el ° 
as a sum of trees in the intermediate field representation. Focusing on the root of the tree, we 
can amputate the unique ciliated vertex into two trivial half-lines (this wont be possible if there 
were wavy loops). Detailing the loop vertex at the other end of the unique wavy line of the tree 
connected to the ciliated vertex leads to the following graphical representation of Y™ el ° 

Melonic 



1PI 





Melonic order 

Figure 12: 1PI melonic two point function 


where we sum over all possible colors for the root wavy line. Combining with (5.1) we get 

<KEc2c') 


^melo t 


^Y Y G ^)-~ 2X Y Y 


(5.2) 


(g c ) 2 + m2 — r™ eZo (g c ) ’ 

^ H c f c f 

where the vector was defined in the previous section. This is a closed equation for the melonic 
self energy. Using Theorem [lj it writes in terms of the function / as 

HT, <:'4) 


Hp) = 2 £ 


gU 2 < c '< 6 


(g 1 ) 2 + m 2 + A^ c , f(q^)' 


(5.3) 


where we recall that = {p,ql,---q g} is a function of p. 

Turning now to the melonic four-point vertex function, we can draw the two end vertices as in 
Figure [13] 



Using the results of section [ 3 ] on the monocolored tree structure of Y™ el ° and taking care of the 
combinatorics we can write 


rr'°(pi,..A) = -4\£[1 - AswWISsmM^, 


(5.4) 
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where M c is define by Tr hc {<j>,<j>) =: E{y l} M p,p 2 ,mp^P >> 




and r yfi?(p) : = E c dint(p c )SymM^ is the simple loop integral with two arcs corresponding to 


Figure 14 



Figure 14: A melonic two point insertion 


Hence 


9int{p ) ^ ^ 


<KEc' <&) 


[(? c ) 2 + m 2 + A Ec' /(?cO 
Using Theorem [l] this means the following relation between / and g 

*(E*£) 


2 • 


9(p) = 4 ( 1 - 53 


(X [t ?) 2 + >n 2 + A Ec- /te) 


(5.5) 


(5.6) 


5.2 Renormalized Equations 


In this subsection we give the renormalized version of the previous equations, 
according to section [4] as 


We rewrite (5.1) 


g 2 (p) = 


<KE c p c ) 


<KE c Pc 


Zp 2 + Z m 2 Tn 2 -r 2 (p) P 2 +m 2 + dzp 2 + 6 m 2 m 2 -T 2 (p) 
Next, we impose the following useful renormalization conditions : 


(5.7) 


A 7 • 1 

' dp 2 lp -° 

(5.8) 

<5m 2 := r 2 (p 2 = 0 ). 

(5.9) 


Note about this expression that, strictly speaking, the function r 2 is a function on Z d , and the 
derivation operation makes no sense. In the last expression, the derivative can be viewed as a new 
function on Z d , obtained from the first by analytic prolongation on (with the preamble remark 
that the function on 7L d admit a natural prolongation on the continuous space), computation of 
the derivative of this new function, and finally restriction to the subset Z d C M d . 
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The renormalized function is therefore obtained from the previous equation and the renormal¬ 
ization conditions, by subtracting its value at p 2 = 0 and its first derivative at the same point: 

rv(p) := r(p 2 ) - r 2 (o) - p 2 ^\ f , s , (s.io) 

A similar equation relates Y^ lo {p) to T™f°(p ), and by the same argument that in the previous 
section, we obtain the following closed equation for the renormalized self-energy: 


r™“>® =-2A r £ £ 


c tf 




c ' t^c 


5 (E d<&) 


_(q c ) 2 + m 2 - T ™?°{q c ) 

d EEdt ') 


(q c ) 2 + m 2 - T ™r lo (qc ) ,Pc ° ^ Pc dp\ ( q c ) 2 + m 2 - T™f°(q c ) lPc ° 

We obtain an equation for f r , the renormalized function / such that 

= -A r £/ r (p c ), 


namely 


f r (.Pc) — 2 Eql,,2<d<6 
KEc' 9c') 


d (Ec' 9c') 


. (9 1 ) 2 + m 2 + A r E^ /r(9c') 
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1 \ \pc—0 Pc . 2/„-l\2 


<KE C ' 9c 1 ') 


dp 2 c (g 1 ) 2 + m 2 + A r E c ' fMl’) 


1 \ lp c =o 


(g 1 ) 2 + m 2 + A r E C ' /r-(9c') 

The renormalized equation corresponding to ( |5.6[ ) follows in the same way. Setting 

= -A ,-£«(!>,=, p')*®) 

c 

(compare with (2.19)) we have the renormalized version of (5.5[) 


Qint, 


,r(P ) = 


C^C 


<KEc'9c') 


<KEc' 9c0 

1 

- \(<?)‘ 2 + m 2 + Ac Ec'/r (9c') 

2 

((f) 2 + m l + K Ec' fr (9c' ) 

2 P=° 


and the renormalized version of (5.6) 


0r(p) = 4 ( 1 - 


<KE C ' 9c') 


^(Ec'9c') 


9 c c Vc L (<f) 2 + m 2 + A r Ec' fMd) (<f ) 2 + m 2 r + X r Ec' /r(9c') 


2 lp=0 


(5.11) 


(5.12) 


(5.13) 


(5.14) 


(5.15) 


(5.16) 


5.3 Existence and Unicity 

The previous closed equations define, in principle, the renormalized mclonic vertex functions. 
Neither the existence nor the unicity of their solutions, however, are obvious at all, since the bare 
equations do not have ultraviolet limit and the renormalized ones typically have zero convergence 
radius in A r because of renormalons (except at very special values such as zero external momenta). 
But we can expand these equations according the multiscale expansion of Section [3] and reshuffle 
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them in terms of the effective amplitudes and effective constants X t of Section |4j Subtractions in 
loop sums such as those of (5.14) and (5.16) will then occur only when the external momentum 


p c has scale strictly lower than the one of q c and the coupling A r will be replaced by the effective 
coupling corresponding to the scale of q c . 

Expanding in a multiseries for all couplings gives therefore an effective expansion with 

• at most ( Ki) n graphs at order n, since as well known, trees proliferate only exponentially in 
their number of vertices, 

• effective melonic amplitudes bounded by ( Kf) n by Theorem [6] (which applies to any effective 
amplitude, hence in particular to the melonic ones), 


effective constants all bounded by the last one A r because of asymptotic freedom (4.37) 


Hence this effective melonic expansion converges and defines a unique solution of the renormal¬ 
ized equations for 0 < A r < (KiK 2 )~ 1 . As usually for flow equations such as (4.37), the solution is 


in fact analytic in A r in a disk tangent to the real axis, with uniform Taylor remainder estimates 
at order s in K s s\ [42]. We leave the details to the reader, but have no doubt that the unique 
solution sum of the effective series is therefore the Borel sum of the renormalized expansion for the 
melonic vertex functions and that this holds not just for N — 1 and 2 but for any number 

2N of external arguments. This completes the control of the melonic sector of the theory: 


Theorem 8 The effective expansions of the renormalized melonic vertex functions converge for 
0 < A r < K _1 to the Borel sum of their renormalized expansions. 

It is tempting to believe that like for tensor models [16], for A r large enough we reach singu¬ 
larities at which phase transitions occur, but this is left to future analytic and numerical study. 


6 Conclusion 

We have studied a simple Abelian TGFT of rank 6 with quartic melonic interactions. We defined 
its intermediate held representation and used it, together with a multi-scale analysis, to prove its 
renormalizability, to compute its beta function and to check its asymptotic freedom. We have 
defined the effective expansion of the model and established uniform exponential upper bounds on 
effective melonic amplitudes. Finally we wrote a closed equation for the melonic approximation to 
the two-point and four-point vertex functions and using the effective expansion we proved that it 
admits a unique solution for small enough stable renormalized coupling. 

Next steps in the analysis of the model might be the numerical analysis of the RG how along 
the lines of H and a full constructive analysis (including the non-melonic sector) of this model. 
The latter would require a non-trivial extension of the techniques of [43], but may be tractable 
thanks to the vector-like nature of the intermediate held. 
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